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Abstract 



The width A^ of the first Bragg's scattering peak in the (111) direction of a 
face-centered-cubic lattice of air spheres can be well approximated by a sim- 
ple formula which only involves the volume averaged e and e 2 over the lattice 
unit cell, e being the (position dependent) dielectric constant of the medium, 
and the effective dielectric constant e e // in the long-wavelength limit approx- 
imated by Maxwell-Garnett's formula. Apparently, our formula describes the 
asymptotic behaviour of the absolute gap width A l for high dielectric con- 
trast 5 exactly. The standard deviation a steadily decreases well below 1% as 
5 increases. For example a < 0.1% for the sphere filling fraction / = 0.2 and 
5 > 20. On the interval 5 E (1, 100), our formula still approximates gap widths 
with a reasonable precision, namely the absolute gap width A^ with a stan- 
dard deviation 3% for low filling fractions up to 6.5% for the close-packed case 
and the relative gap width A r L from 4.2% to 8%. Differences between the case 
of air spheres in a dielectric and dielectric spheres in air are briefly discussed. 



(3. Opt. A: Pure Appl. Opt. 1, 471-475 (1999); 
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I. INTRODUCTION 



The propagation of light in a periodic dielectric medium has recently attracted much 
attention due to the possibility of opening a gap in the spectrum of electromagnetic waves 
for both polarizations and all directions of the incident waves |l|-|5||. In such a medium, 
the density of states (DOS) can be, in a certain frequency interval, either reduced down 
to zero (photonic band gap) or enhanced with respect to the vacuum case. The changes 
in the DOS affect various physical quantities. The most transparent is the change in the 
spontaneous emission rate of embedded atoms and molecules which may have applications 
for semiconductor lasers, hetero junction bipolar transistors, and thresholdless lasers [JTJ] or 
to create new sources of light for ultra-fast optical communication systems. 

Existence of the full photonic band gap was first demonstrated at microwaves ||. Re- 
cently, thanks to the intense experimental effort, we have witnessed a significant progress 
in fabrication of complete photonic-bandgap structures at near-infrared In two di- 

mensions, complete photonic-bandgap structures have been fabricating for TM polarization 
||. In three dimensions, a promissing structure has been fabricated by Sandia's group J7j. 
Nevertheless, in one direction this structure extends less than two unit cells and there is 
ongoing experimental search to improve its properties. One of the most promising candi- 
dates to achieve a complete photonic bandgap at optical wavelengths and fabricate large 
enough structures at near-infrared wavelengths are collodial systems of microspheres. In- 
deed, the latter can self-assemble into three-dimensional crystals with excellent long-range 
periodicity with the lattice constant well below infrared scale PHTBl. This long-range peri- 
odicity gives rise to strong optical Bragg's scattering clearly visible by the naked eye and 



already described in 1963 [14]. Monodisperse collodial suspensions of microspheres crystal- 



ize either in a face-centered-cubic (fee) [|TTJ[T5|[T6| or (for small sphere filling fraction) in a 
body-centered-cubic (bec) lattice |12|]. Using suspensions of microspheres of different sizes 
one can also prepare crystals with a complex unit cell (containing more than one scatterer). 
Both the case of "dense spheres" |T9 and "air spheres" |I7j when the dielectric constant of 
spheres e s is greater and smaller than the dielectric constant £& of the background medium, 
respectively, can be realized experimentally. There is a significant difference between the two 
cases, since, according to numerical calculations, simple dielectric lattices of homogeneous 
spheres [P~5| — PC|] in air do not exhibit a full photonic band gap, while for air spheres a full 
band gap can open for a simple fee lattice [|I8| -|2"0"11 ■ Unfortunately, the required dielectric 

contrast 5 = max(e s /e b , Eb/e s ) for opening the full band gap, either ~ 8.4 obtained using 

the plane wave method |19|], or, ~ 8.13 obtained by the photonic analoque of the Korringa- 



Kohn-Rostocker (KKR) method ]^U,^T|, is currently out of experimental reach at optical 
and near-infrared frequencies for photonic colloidal structures. The absence of a full gap in 
this frequency range in currently available collodial crystals of homogeneous and single size 
spheres does not mean the absence of interesting physics in this weak photonic region. For 
example, the change in the spontaneous emission rate of dye molecules in an fee collodial 
crystal can be observed already at a relatively low 5 ~ 1.2 p2| . 

In contrast to the full gap, Bragg's reflection can be observed for arbitrarily small 5 as 
long as a sample has sufficient long-range periodicity. Analysis of Bragg's scattering might 
not only be useful to understand the physics of photonic crystals, but it has already found 
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practical application in distributed feedback lasers in the visible region of the spectrum |23 . 
The first Bragg's peak can be characterized by the width of the (lowest) stop gap (gap at a 
fixed direction of the incident light) at a certain point on the boundary of the Brillouin zone. 
We focus here on the case of a simple fee lattice of air spheres |17J, which is among the most 



promising candidates to achieve a full photonic band gap. For an fee lattice, it is convenient 
to consider Bragg's scattering in the (111) direction which corresponds to the L direction 
of the Brillouin zone (see for the classification of special points of three-dimensional 



lattices). Apart from numerous experimental data now available 151,16, 2? |, there are at 
least two other reasons for this choice. First, the width of the first stop gap takes on its 
maximum at the L point and, second, experimental techniques make it possible to allow one 
to grow collodial crystals such that the L direction corresponds to normal incidence on the 
crystal surface. 

Let e(r) be the dielectric constant of an fee photonic crystal. One has e(r) = e s if r is 
inside the sphere and e(r) = e b otherwise. Let / be the sphere filling fraction, i.e., volume 
of the sphere(s) in the unit cell per unit cell volume. Once / is fixed, the spectrum is only 
a function of the dielectric contrast 5. By a suitable rescaling, one can always set e s — 1 for 
the case of "air" spheres (e b = 1 for the case of "dense" spheres). As 5 and / are varied, 
both the absolute L-gap width Ax and the L-midgap frequency v c change. As a function 
of S, Ax, (5) takes on its maximum at some 5 = S m (f) while v c (5) monotonically decreases. 
We address the question of whether the width Ax can be understood in terms of simple 
quantities, namely, the volume averaged dielectric constant, 

£ = fe s + (1 - f)e b , 

the volume averaged £ 2 {r), 

and the effective dielectric constant e e ff. The latter characterizes optical properties of the 
crystal in the long-wavelength limit and is (theoretically) determined by the slope of the 

— 1/2 

linear part of the band structure, e e J = lim^o duj/(cdk). Note that due to the vector 
character of electromagnetic waves, e e ff differs from e, in contrast to the scalar case where 



e e ff = e p4j| . One can show that for any /, e b , and e s , 

e ef f <e<^. (1) 

Equality in ([!]) occurs if and only if either / = or / = 1, or, if 8 = 1, i.e., if e b = e s . The 
effective dielectric constant can be well approximated f2C||24|l by Maxwell-Garnett's formula 



eeff*e%f G = e b (l + 2fa)/(l-fa), (2) 

where, for a homogeneous sphere, the polarizability factor a = (s s — e b )/{e s + 2e b ). Note, 
however, that in the case of air (dense) spheres slightly overestimates (underestimates) 
the exact value of e e ff as calculated from the band structure [^0| , ^| . 
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II. RESULTS 



Obtaining exact analytic results for dielectric lattices turns out to be notoriously difficult 
and numerics has been the main tool to understand photonic gaps so far [^,^|,p!8 H20|j24| . A 
simple analytical formula, if any, may be a good starting point for obtaining a better insight 
into the problem. It was rather surprising to find out that such a formula can be found for 
the L-gap width A^. Namely, in the case for air spheres, A^ can be approximated by the 
formula (see Fig. [I]) 



A L ^Cg=C 



■eff 



1/2 



£. 



(3) 



For a given filling fraction /, the constant C = C(f) was determined by taking the average 
over Ai/g where A^ is the L-gap width calculated numerically using a photonic analogue 
[21 1 of the KKR method |28 [. The latter method gives results which are in excellent agree- 



ment with experimental values [16| . 

Apparently, for sufficiently high 5 ^> S m (f), our formula captures the asymptotic be- 
haviour of the absolute gap width A l exactly. The standard deviation a steadily decreases 
well below 1% as one investigates region S > S c for higher and higher S c . For 6 > 20 one has 
a < 0.1% for / = 0.2, a < 0.4% for / = 0.1, and a < 0.3% for / = 0.4 (see, for example, 
Fig. Q). For 5 > 36 and / = 0.6, a < 0.5%, while for the close-packed case, a < 1%. If, 
however, in the latter case 5 > 50, a drops below 0.7%. 

For 5 G (1, 100], our formula (p3) still describes A^ with a reasonable accuracy ranging 
from 3.3% to 6.5% (depending on the filling fraction). The values of C, their standard 
quadratic deviation a, and the relative error a r = a/C are shown in Tab. I. Approximately 
thirty values of the dielectric contrast within the interval 5 G (1, 100] were taken for every 
filling fraction considered. For a given filling fraction, the main part of the error is picked 
up 

TABLE I. The values of C, their standard quadratic deviation a, and the relative error 
cr r = a/C for different filling fractions and 5 G (1, 100]. 





/ = 0.1 


/ = 0.2 


/ = 0.4 


/ = 0.6 


/ = 0.74 


c 


0.762 


0.868 


0.875 


0.808 


0.736 


a 


0.031 


0.03 


0.029 


0.038 


0.048 


o r 


4.1% 


3.4% 


3.3% 


4.7% 


6.5% 



around 5 = S m (f) for which A^ takes on its maximum. At the maximum of the L-gap width 
A i, our formula (0) gives persistently a slightly lower value for A L . Note that for moderate 
5 Maxwell-Garnett overestimates e e ff for the case of air spheres [2"0" , 2"1 |. Therefore, using the 
exact e e ff may reduce errors further. 

According to Tab. I, the quantity C shows a weak dependence on / which can be 
approximated with high accuracy (relative error 2.5%) by the formula 

C(f)=C + 0.Uf(2f m -f)/ti 



c2 
1 m" 



(4) 
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Here Co ~ 0.74 is the minimal value of C and f m is the filling fraction for which C takes on 
its maximal value. Table I indicates that C takes on its minimal value Cq at the extreme 
filling fractions / = and / = 0.74, and its maximal value is C m ~ 0.88 at f m ~ 0.74/2. 
The factor 0.14 in the interpolation formula (|J) is the difference C m — C . Using C(f) in @ 
does not raise the relative error a r more than 0.4% for intermediate filling fractions. Fig. 
|I] shows approximations to the L-gap width for / = 0.4 using formula (§) with optimized 
C taken from Tab. I and with C(f) given by the formula (ffl). As the dielectric contrast 5 
increases, first increases to its maximal value and then slowly decreases as 5~ 1 ^ 2 . This 
behavior is well reflected by our formula (^) which in the limit 5 ^> 1 yields 

Al ~ C(f) . 2 (tt^== " lV /2 ^ 1/2 - (5) 

/(l-/)(2 + /) V2v/W ; 



Since the L-midgap frequency v c changes as / and S are varied, an invariant characteristic 
of Bragg's scattering is provided by the relative L-gap width 

A r L = A L /v c . 

A r L increases monotonically as 6 increases (see Fig. ^) and saturates very fast for 5 > 5 m . 
Our observation here is that the L-midgap frequency v c can be well approximated by 



u c ^ck L /(2rcn^ f G ), (6) 
where n^P = \ s^P and ki, is the length of the Bloch vector at the L point. In units where 



the length of the side of the conventional unit cell of the cubic lattice [29] is A = 2, one has 



hi/it = V0.75. Recent measurements of v c for moderate 5 \TG] agree well with formula (||) 
(see also Fig. |2]). For all filling fractions considered 5 within the range and 1 < 5 < 100, 
the maximal deviation of the L-midgap frequency given by formula @ is less than 8% with 
respect to the exact value. Therefore, the formula 

Al « 2^A L /k L (7) 

is a natural candidate to describe A r L . However, as shown in Fig. ^ formula (jg) systemati- 
cally overestimates the exact value of v c by a little bit. This systematic error is also apparent 
from Fig. |j. Due to the systematic error, the relative L-gap width is described by the 
formula (|7|) with a slightly larger relative error than is Al by the formula (0). There are 
now two main contributions to the errors, one around the maximum of A^ and the other 
due to the systematic error. However, in the asymptotic region 5 ^> 5 m the first contribution 
disappears whereas the systematic error saturates (see Fig. 0). As a result, in the asymptotic 
region, the relative error a r is still within w 1%. Fig. § shows that even at 5 = 100 the error 
in A ^ is less than 2%. 

For 5 G (1, 100], our formula (f|) still describes A^ with a reasonable accuracy ranging 
from 4.1% to 8% (depending on the filling fraction) For a given filling fraction, the relative 
error a r = a/R was determined by calculating the standard quadratic deviation a of the 
average value R of the ratio R = /\ r £ exact j / ^'n approx } where /\ r £ exact j s the exact value of A^ 
calculated numerically and /\ r £ a PP rox j s n s approximation calculated using Eq. (0). The 
values of R, a, and a r are collected in Table II. 
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TABLE II. The average value R, the standard quadratic deviation a, and the relative error 
a r for A r L approximated by Eq. ([71) for different filling fractions and 5 £ (1, 100]. 





f = 0.1 


/ = 0.2 


/ = 0.4 


/ = 0.6 


/ = 0.74 


R 


1.004 


1.012 


1.025 


1.030 


1.031 


o 


0.044 


0.042 


0.052 


0.068 


0.082 


a r 


4.37% 


4.17% 


5.05% 


6.61% 


7.95% 



One expects a deviations [of the order 5% from the behavior described by the formula (^)] 
only in a rare case when a Mie resonance crosses the edge of the L-gap |IIJ . 



III. DISCUSSION 



Formulas (|) [together with (g)], (1), and (0) are the main results of this work. They fit 
nicely experimental data on Bragg's scattering in fee photonic crystals of air spheres p~5^[T7[| . 



Note that also characterizes the transmission of light through such a crystal (see |27| for 
microwaves). The fact that such simple relations can describe one of the photonic gaps has 
been completely unexpected. Indeed, the numerical calculation of photonic band structures 
is a great deal more involved than that in the case of scalar waves (including the case of 
electrons) where no analog of formulas @ and (^) is known. Numerics has been the main tool 
to understand photonic gaps |]3|J4]JT8|- pO| , p4l . This is also the case of two recent discussions 
of Bragg's scattering in the (111) direction ||32|| . A previous attempt to understand Bragg's 



scattering in photonic crystals involved an introduction of a "photonic strength" parameter 
^ = 3fa \W\. It was shown that the dynamical diffraction theory ||33|| , which is well known 
in x-ray diffraction, already fails to describe Bragg's scattering in a photonic crystal for 
* « 0.5 |T§. 

Formulas (^) and (^) immediately raise questions whether one can understand and derive 
them analytically. The L-gap width for fee structures is a natural measure to characterize 
their scattering strength, because, in contrast to the full band gap, ^ for arbitrarily 
small / and 5. The latter suggests to take A r L given by Eq. (0) as a natural "photonic 



strength" parameter for the air spheres case. Neither the parameter \& [16|, nor the parame- 



( o _o\ 1/2 



ter e r = (e 2 — e 2 J ' j e, introduced in [18], are directly related to a gap width. However, it 

turns out that formulas (|) and (|7|) cannot be applied to the case of dense spheres. The sim- 
ple fee lattices of air and dense spheres have for the same dielectric contrast rather different 
behavior with respect to the full photonic band gap [TJ-F2U] and to the first Bragg's peak 
|n| , ^5H l7| . Our numerical calculation shows that, for dense spheres, A r L does not increase 
monotonically with 5 as in the case of air spheres. Instead A r L first reaches a local maximum, 
then it returns to zero and only afterwards starts to increase monotonically |31| . This behav- 
ior is reminiscent to that of the relative X-gap width (X is another special point of the 
Brillouin zone of an fee lattice [53]) in the case of air spheres Pj2^|. It has been argued that 



the vanishing of A^ is due to the vanishing of the scattering form-factors P,p7|. Also, if A r L 
is plotted against the filling fraction, one observes that the maximum of A^ shifts to lower 



/ for dense spheres and towards close-packing for air spheres [[Hp. It would be interesting 
to understand what causes this different behavior. The latter can be partially attributed 
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to the fact that, for a lattice of spheres, e& no longer describes the dielectric constant of 
the surrounding medium, which is instead described by the effective dielectric constant e e ff. 
Therefore, the bare dielectric contrast 8 is renormalized to 8 e ff = max(£ s / £ e ff, £ eff/ £ s), where 
1 < e e ff < 5 for e s 7^ £&. Given the bare dielectric contrast 5, one finds that the renormal- 
ized dielectric contrast = e s /e e ff m the case of dense spheres is always smaller than 
the renormalized dielectric contrast 5^* = e e ff in the case of air spheres ||20|| . The latter is 
easy to verify in the limit when the bare dielectric contrast 5 tends to infinity, where the 
Maxwell-Garnett equation (|2]) implies 

5 d eff -5(1- /)/(l + 2f) < 5 a eff -5(1- f)/(l + f/2). (8) 

Nevertheless, a full understanding of the differences between the lattices of air and dense 
spheres still remains a theoretical challenge. 

IV. CONCLUSION 

To conclude we have found that, despite of the complexity of the problem of propagation 
of electromagnetic waves in a periodic dielectric medium, the absolute and the relative 
width of the first Bragg's peak in the (111) direction for an fee lattice of air spheres can be 
accurately described by the simple empirical formulas (H) and (^), respectively. Apparently, 
for sufficiently high 5 > 1, our formula (|3|) captures the asymptotic behaviour of exactly. 
Indeed, the relative error a r steadily decreases as one investigates region 5 > 5 C for higher 
and higher 5 C . For all filling fractions a r falls well below 1% if sufficciently high 5 is taken. 
For example for 5 > 20 one obtains a r < 0.1% for / = 0.2 and a r < 0.3% for / = 0.4. 
For 5 G (1, 100] formula (|3|) still describes with a reasonable precision, namely, with 
the relative error ranging from 3.3% to 6.5% (depending on the filling fraction). The main 
contribution to the error is picked up around 5 = 5 m (f) for which A^ takes on its maximum. 
At 5 = 5 m , our formula (0) gives persistently a slightly lower value for A L . The relative 
L-gap width A^ is described by the formula (0) with a slightly larger relative error ranging 
from 4.1% to 8% (depending on the filling fraction). The reason is that there are now two 
main contributions to the error, that around the maximum of A^ and the second systematic 
error due to the overestimation of the L-midgap frequency v c when using Eq. ([]). All the 
formulas only involve the effective dielectric constant of the medium e e ff approximated by 
Maxwell-Garnett's formula (0), and volume averaged e(r) and £ 2 (r) over the lattice unit cell. 
Since e, e 2 , and e e ff have well-defined meaning for any lattice, this suggests that a similar 
gap behavior may occur for other lattices. It would be interesting to find out if the same is 
true for the width of the full photonic band gap |yj • 
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FIGURES 
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Refractive index contrast n b /n s 

FIG. 1. The L-gap width in units c/ttA, where c is the speed of light in vacuum and A is a 
length of the side of the conventional unit cell of the cubic lattice [25], for an fee lattice of air 
spheres with / = 0.4 calculated exactly (the solid line) and approximated by the formula (||) with 
optimized C (the dashed line) and with C given by the formula (Q) (the long-dashed line). The 
last two curves almost overlap. 
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Air spheres 




- 1 o.92 ' ■ 1 ■ 1 ■ 1 ■ 1 

1.0 3.0 5.0 7.0 9.0 

Refractive index contrast n b /n s 

FIG. 2. The L-midgap frequency v c for air spheres in the units fei/(27rn^ G ), where ki is the 
length of the Bloch vector at the L point and n^j^ is the effective refractive index of the medium 
calculated by Maxwell-Garnett's formula (||). The dot-dashed line is for the sphere filling fraction 
/ = 0.1, the long-dashed line is for / = 0.2, the dashed line is for / = 0.4, finely dotted line is for 
/ = 0.6, and the solid line corresponds to the close-packed case (/ = 0.74). 
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Air spheres 



t ■ 1 ' 1 1 r 




Refractive index contrast n b /n s 

FIG. 3. The relative L-gap width (the L-gap width divided by the midgap frequency) for an fee 
lattice of air spheres shows a rather simple dependence on the refractive index contrast rib/n s = y/S. 
The dot-dashed line is for the sphere filling fraction / = 0.1, the long-dashed line is for / = 0.2, 
the dashed line is for / = 0.4, finely dotted line is for / = 0.6, and the solid line corresponds to 
the close-packed case (/ = 0.74). 
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Air spheres 
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FIG. 4. The relative L-gap width A r L for the close-packed fee lattice of air spheres (solid line) 
and its approximation using the formula (^) with optimized R (long-dashed line). Once R is 
optimized, a r can be reduced to ~ 1%. Even at 8 = 100, the error in A r L is less than 2%. 
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